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Numerical Solutions of Supersonic and Hypersonic Laminar 

Compression Corner Flows 

C.M. Hung* and R.W. MacCormackf 
NASA Ames Research Center, Moffett Field, Calif. 


An efficient lime-splitting, second-order accurate, numerical scheme is used to solve the complete Navier- 
Stokes equations for supersonic and hypersonic laminar flow over a two-dimensional compression comer, A 
fine, exponentially stretched mesh spacing is used in the region near the wall for resolving the viscous layer. 
Good agreement is obtained between the present computed results and experimental measurement for a Mach 
number of 14.1 and a Reynolds number of 1.04x 10* with wedge angles of 15% 18% and 24% The details of the 
pressure variation across the boundary layer are given, and a correlation between the leading edge shock and the 
peaks in surface pressure and heat transfer is observed. 


1. Introduction 

C ONTINUING advances in numerical methods and com- 
puter capabilities have now made feasible many flowfield 
calculations which were formerly intractable. One such 
problem which has renewed interest is supersonic or hyper- 
sonic flow over a two-dimensional compression corner. This 
problem has received considerable attention within the past 
decade because of its importance to the design engineer in 
predicting the pressure and heat loads at a wing-flap junction 
on re-entry vehicles. When flow separation occurs, reduced 
flap effectiveness results, and in some regions, the surface 
heating could become severe on a maneuverable re-entry 
vehicle. 

The problem to be considered is illustrated schematically in 
Fig. I. The pressure rise generated by the wedge extends up- 
stream along the flat plate, thickening the boundary layer, 
and results in a complicated interaction between the viscous 
flow near the body surface and the outer inviscid stream. Sin- 
ce the inner part of the boundary layer may not have suf- 
ficient momentum to overcome the combined effects of skin 
friction and adverse pressure gradient, the interaction can 
lead to flow separation for certain ranges of Mach number, 
Reynolds number, and wedge angle. The separated boundary 
layer will then become a free shear layer external to a recir- 
culating inner flow near the corner. Reattachment occurs 
because of the interaction between free shear flow and the 
outer flow. The surface pressure continues to rise through the 
separated and reattached regions, until the boundary layer 
reaches a minimum thickness or ”neck.” Downstream of the 
neck, the boundary layer returns to a normal state of weak in- 
teraction with the outer inviscid stream at a new Mach num- 
ber. Although, in most practical situations the region of shock 
wave and boundary-layer interaction is turbulent, at high 
altitude flight, fully laminar flows can exist and are important 
for design considerations. 

Previous theoretical treatment of such a problem has 
usually been made with the boundary-layer equations together 
with a “coupling” equation relating the development of the 
inner viscous flow to the outer flow. The governing partial 
differential equations can then be solved by finite difference 
techniques*'^ or can be expressed as integral relations and 
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solved as ordinary differential equations. ^ This treatment, in 
general, involves the question of uniqueness because certain 
portions of the flowfield contain substantial upstream in- 
fluence, and the initial and downstream boundary conditions 
cannot be completely specified. In some integral tech- 
niques,'*'^ there is also the question of so-called jump con- 
ditions for supercritical to subcritical types of boundary 
layers. 

Consideration of the Navier-Stokes equations avoids these 
questions and removes some restrictive assumptions, viz., that 
the static pressure is constant across the boundary layer, and 
that the viscous and inviscid flows interact only along a line at 
or near the edge of the boundary layer, which can be difficult 
to define in hypersonic flow. Also, as mentioned by Van 
Dyke,^ a solution of the Navier-Stokes equations is necessary 
in the immediate vicinity of a sharp corner. 

Carter^ has obtained numerical solutions of the Navier- 
Stokes equations for laminar flow past a compression corner 
at low Mach numbers. Fle used the Brailovskaya difference 
scheme, which is first-order accurate in time and second-order 
accurate in space. In the present study a more efficient 
numerical method, which is second-order accurate in both 
time and space, is used to calculate supersonic and hypersonic 
cases. 

The complete time-dependent Navier-Stokes equations are 
split into two sets, one for the x-derivative and the other for 
the j-derivative. The advantage of the split system is that the 
computation proceeds with larger time increments because the 
stability criterion is less stringent. Each calculation starts 
from a uniform flow with appropriate boundary conditions. 
The boundary layer and shock waves generated by the leading 
edge and compression corner develop with time, until a steady 
state is achieved. A fine, exponentially stretched mesh spacing 
is employed in the region near the wall for resolving the 
viscous layer, and a coarse uniform mesh spacing is used in 
the outer region where viscous effects are negligible. The 
solution is advanced more frequently in the fine mesh region, 
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Fig. 1 Schematic diagram of supersonic flow over a sharp flat plate- 
wedge compression corner. 
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using smaller time steps, than in the coarse mesh region. Com- 
putational efficiency is further enhanced by using elongated 
cells (nonisotropic mesh spacing). 

A novel, recently devised fourth-order damping term is in- 
corporated to maintain numerical stability in the presence of 
steep pressure gradients, and hence, fine mesh spacing near 
the leading edge, required by previous studies, is avoided. 
Moreover, high Mach number cases with strong induced 
shocks may be solved. The results are in good agreement with 
Carter’s computed data, with significantly less computation 
time for a wedge angle of 10°, Mach number of 3.0, and 
Reynolds number of 1.68x10^, based on freestream con- 
ditions and the distance from the leading edge to the corner. 
For a severe test, the experiments conducted by Holden and 
Moselle^ are selected for comparison at Mach number 14H 
and Reynolds number 1.04 x 10^ with wedge angles of 15 , 
18°, and 24°. For the large-wedge-angle case, the pressure 
rises by a factor of 50 across the interaction. Good agreement 
is obtained between the present computed results and ex- 
perimental measurements of wall pressure, skin friction, and 
heat transfer. 

With the whole flowfield calculated, we are able to study 
the details of the shock wave and viscous flow interaction. 
The most striking feature is the large, normal static pressure 
gradients within the boundary layer, which contrasts with the 
classical boundary-layer assumption. Instead of starting with 
a local similarity flat plate solution, the generation of a 
leading edge shock provides another interesting feature. For 
the experimental cases of Holden and Moselle^ we observed 
from the calculated flowfields that the interference of the 
leading edge shock with the induced shock significantly af- 
fected the pressure and heat transfer distributions. 

II. Analysis 


Governing Fqualions 

The time-dependent Navier-Stokes equations, in two 
dimensions, neglecting body forces and heat sources, may be 
written in integral form as 

— { Ud\o\+[ HndS = 0 (1) 

dt J volume surlacc 
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the volume of integration. The equation of slate relates the 
pressure p and density p to temperature T and specific internal 
energy e. The perfect gas relations arep = p/?T and e = C,T, 
The viscosity coefficient p is assumed a function of tem- 
perature only, and is evaluated by Sutherland’s semiempirical 
formula p = 2.270 x 10 T^ “/( T-i- 198.6) (slug/ft-sec). For 
the conditions of the experiment and the calculated tem- 
perature distribution we assume that the above relations are 
adequate for this study. The coefficient of thermal con- 
ductivity K is computed by assuming a constant Prandtl num- 
ber Pr and specific heat c ^ , so that K = yC* p/Pr. 


Time-Split Numerical Technique 

A time-split, two-step MacCormack scheme is used which 
solves the set of two-dimensional equations by reducing them 
to two sets of one-dimensional equations. The system of 
equations can be solved in one orthogonal x, y coordinate 
system for an arbitrary quadrilateral volume element, as 
described by Deiwert. The typical element used is shown in 
Fig. 2, where B is the inclination of the parallelogram element 
with respect to the x-coordinate. The difference equations ap- 
proximating Eq. (1) are then written as follows 
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and e,,6y are unit vectors of the orthogonal x,y coordinate 
systern, and h a unit normal vector of the surface enclosing 


The subscripts /, j refer to a special mesh of point with 
spacing Av, and Ayj, and the superscript n refers to times 
i = nAt where At is tl^ tim£_increment that the solution is ad- 
vanced. The bar on Fand G indicates that the predicted quan- 
tities 17 are to be used in the evaluation of these terms. Note 
that the operator Lx accounts for the convection and stress 
terms acting at the vertical sides of the typical mesh element 
show n in Fig. 2 and the Ly operator accounts for the terms ac- 
ting on the inclined sides. The dot product of Eq. (1) is FAyj 
for a vertical side and (G + Ftan 0,)Axi for an inclined side. 
The Lx and Ly operators are applied so that Eq. (1) is ap- 
proximated at each time-step for each element of the non- 
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Fig. 2 Parallelogram volume element used in computation. 
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Fig. 3 Computational boundaries (above), and mesh point and cell 
system (below). 


orthogonal computational mesh. Although the operator 
sequence of Ly (AO Lx (AO is accurate only to the first order, 
MacCormack has shown that symmetrical operator sequences 
such as 

'j 

are accurate to the second order both in time and space. 

To evaluate the viscous and heat conduction derivatives for 
the nonorthogonal cell (Fig. 2), the differences used are as 
follows: Using 0 as a dummy variable, for the Lx operator 
predictor 

\dx hj3 x,^,-x, 

\dy ) ij yj+,-yj-, 

and for theL>^ operator predictor 

V 9X / IJ ^i+i —Xj^i yj+i~yj 

/ d4> \ = 

where 

Aij=Vi (tan 0,^./-i-tan 6^) 

Bij^V4 (tan tan 0, + tan 

Similar approximations are used in the corrector steps, such 
that the overall treatment of derivatives results in centered 
second-order-accurate approximations. For variable mesh 


Reynolds number is computed by {^y),,,u^ = (Lv7?Cy ). The 
leading edge is placed at x = 0 and the distance from the 

leading edge to the corner, is used as the characteristic length. 

The appropriate boundary conditions are shown in Fig. 3. 
At the upstream boundary a few mesh points ahead of the 
leading edge, the flow is assumed uniform at supersonic 
freestream conditions ^) . The downstream 

boundary is positioned far enough from the corner so that all 
the gradients in the flow direction can be set to zero. Though 
this condition is not exact, the boundary layer in the vicinity 
of this exit is parabolic, and the remainder of the flow is 
supersonic; hence, it is not expected that this condition will in- 
troduce significant error in the region of interest upstream. 
The upper boundary is specified by the freestream conditions 
ahead of the leading edge shock and by simple wave ex- 
trapolation downstream of the shock. The wall surface is 
assumed impermeable, and nonslip boundary conditions are 
applied. The wall is treated either as isothermal or adiabatic, 
and the wall pressure is evaluated from an approximation of 
the j^-momentum equation at the wall. During the calculation 
of the inner mesh, transport and stress at the internal boun- 
dary are saved, and their average net quantities are then 
used as boundary conditions for the outer mesh flowfield. 

Computational Time-Step 

The maximum time-step for which the calculation will be 
stable is determined by the CFL and viscous stability 
requirement. By applying a von Neumann stability analysis 
separately to the linearized inviscid and viscous parts of Eqs. 
(1)'’ we arrive at the following criteria: For the Lx{At ) 
operator ' 

\u\ +ff+ — ( - ) ~|] 

L p Ax Ay, JJ . 

J min 

and for the Ly ( A/^, ) operator 




spacing the resulting difference equations are second-order 
accurate in the computational coordinate system. 

Mesh and Boundary Conditions 

Figure 3 shows the computational domain and the mesh 
point and cell system used for calculation. Note that the in- 
clination angle B is zero on the flat plate and is constant along 
the wedge. The mesh is equally spaced in the jr-direction, but 
they-direction, a fine, exponentially stretched mesh spacing is 
used in the region 0.0<y<hj near the wall for resolving the 
viscous layer, and a coarse, equally spaced mesh is used in the 
outer region hj<y<h, where viscous effects are negligible. 
The minimum spacing for the fine mesh is dependent upon the 


I Ay^ / \v~~u tan 0,\ -ya sqcO,-\- — ( -I- — ^ ) \ 

^ ^ P Ay. Ax JJ . 

^ mm 

where a is the local sound speed. These stable time-steps are 
larger than those of explicit finite difference methods without 
time splitting. Four different time-steps are used in 
calculations: two for the Lx operator A/^.,., At^j and two for 
the Ly operator Ai^j^ in the coarse and fine mesh, respec- 
tively. In general, for the present calculations At^j >At,j and 
Afy^.>Af ^^.. Hence, the computational efficiency is enhanced 
by use of the operator sequences for the fine mesh region 

ILyiAtj- /2)Lx{Aty)Ly{A(y /2) ] 
and for the coarse mesh region 

lLx{At,)Ly {2 At,) Lx (A/,)] 
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where 

A/y = min 2Afyf] 

At^ = m‘Atf<min[A(j,^~, ViAty^.] 

and m is an integer, representing the number of applications 
of the operator sequence in the inner mesh region for each 
time-step which the solution is advanced in the outer 
region. As m becomes large, most of the computing time is 
spent in the calculation of the fine mesh-region. 

Nonlinear Instability 

Several types of nonlinear instability are encountered in the 
present calculation. The remedies for these have been 
previously discussed by MacCormack,*^ Baldwin and Mac- 
Cormack,^-* and Kutler, Sakell, and Aiello. ''' One of them oc- 
curs, for example, in the convective flux across a mesh surface 
w hen there is an expansion in which the velocity normal to the 
mesh surface changes sign. When this condition occurs, in- 
stability can be avoided by employing the average normal 
velocity. Another nonlinear instability can be removed by ad- 
ding a product fourth-order damping term AUij to the right- 
hand side of the predictor and corrector equation of each 
operator. For the Lx operator 



Fig. 4 Comparison of wall pressure and skin friction for =3.0, 
Re; =1,68x10^, a =10". 

pressure, skin friction, and heat transfer along the wall. The 
flow conditions selected are as follows: 

Rei^ = 1.04x]0^ Pr^O.72 

T^ = J30°R L=^}.44h T^ = 535''R 

The wall is highly cooled, so that the ratio of wall temperature 
to the adiabatic wall temperature is 0.1 18. Note that. 
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where Cj, and Cy are two assigned constants with values bet- 
ween 0.0 and 0.5. The damping form of the damping term is 
more compact than standard fourth-order smoothing terms. 
It requires data at just three mesh points instead of the usual 
five points. For shock wave calculations this compactness 
achieves better shock resolution. In the present investigation 
these terms played important roles in smoothing the ripples in 
the solution in the region near the leading edge and corner 
shocks. With this smoothing treatment, fine-mesh spacing 
near the leading edge, such as that found necessary by 
Carter^ as well as others, in the x-direction is avoided. This 
reduces considerably the computation time of the present 
method from that required by the other investigators. Even 
more important, high Mach number cases with strong shocks 
can be solved without difficulty. 

III. Results and Comparisons 

The first case computed was for =3.0 supersonic How 
over a 10° wedge with an adiabatic wall. The flow conditions 
correspond identically to one of the cases studied by Carter,^ 
and are given as follows: 

M^=3.0 Re;, = 1.68 X 10^ a=/0° 

T^=390R L=2.4 in. T^, = T^U {y~ I /2)Ml^ 

Here, wall temperature, is assumed equal to the 

freestream stagnation temperature. For this case 86 mesh 
points were equally spaced in the x-direction, and 28 mesh 
points were used in the j-direction, 15 for the inner mesh and 
13 for the outer mesh. Figure 4 shows the present computed 
surface pressure and skin friction distribution. The skin fric- 
tion coefficient is defined as Cy= (r„,/poowi)* These 
distributions are in good agreement with the results of Carter. 

The experiments selected for comparison were conducted 
by Holden and Moselle.® Measurements were made of 


in this case, both the Mach number and Reynolds number are 
almost one order of magnitude higher than for the previous 
calculation. The computational domain extends from (x/L) 
= -0.077 to (x/L)= 1.859, with the fine-mesh boundary 
placed at hr = 0.096 ft and the outer boundary at h = 0.2S ft. 
In the x-direction 90 mesh points were used with a mesh 
spacing of Ax = 3.130 x 10 ft. In the >'-direction 30 mesh 
points 20 for the fine mesh and 10 for the coarse mesh, were 
used. Mesh spacing varied from Ay = 2.087 x 10“^ ft to Ay 
= 9.842x 10“^ ft for the fine mesh; for the coarse mesh, 
equal spacing of Ay = 2.024 x 10"^ ft was used. In general, the 
stable time-steps were about A/ 3 ,y^ = 3.297 x 10“^ sec, A/^ 
= 3.216x10-® sec, A/^^ = 3.272 x 10 ‘® sec, A/y^ = 1.300x 
10"^ sec, and hence, m = \0. The freestream velocity was 
=7.880 X 10^ fps, and it took about 200 time-steps (or 1.3 
msec) to reach a steady state. During this time the freestream 
traveled a distance seven times the length of the flat plate. 
About 32 min on a CDC 7600 were required to perform the 
calculations. 

The calculation was first made for a plate plate. Shown for 
comparison in Fig. 5 is the wall pressure distribution of the 
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Fig. 5 Comparison of 
strong interaction 
theorj^ with present 
results for wall pressure 
distribution on flat 
plate. 
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Fig. 6 Comparisons of: a) wall pressure and b) heat transfer and skin 
friction along wall for cr = IS"*. 


••I X,^g||/L 

Fig. 8 Comparisons of: a) pressure and b) heat transfer and skin fric- 
tion for a = 24*’. 




b) 


Fig. 7 Comparisons of: a) pressure and b) heat transfer and skin fric- 
tion along wall for a = 18®, 


present computed result and the result theoretically predicted 
by the strong-interaction analysis of Bertram and Black- 
stock.^^ The analysis was based on the hypersonic similarity 
theory that 

p/Poo = 0.8J-hy4(y(y^/)/2}'^^ ffxoo 

where 


1,648 


y-1 


-\- 0 . 352 \y Xoo = 


VRe, 


Here, Xcx is a hypersonic parameter; /3 is a parameter related 
to the wall temperature; and C, the Chapman-Rubesin con- 
stant, is evaluated as C=T^^(T^.) /T^^{T^) , Excellent 
agreement is obtained between the theoretical prediction and 
the numerical result. Also shown in Fig. 5 is an experimental 
value of Holden and Moselle.^ The reason for the 
disagreement between the present result and the experiment is 
unknown. Nevertheless, as will be seen later, good agreement 
was observed downstream in the interaction region. 

Figures 6-8 show the detailed comparisons of the present 
computed results with the experimental measurements for 
pressure, heat transfer, and skin friction for three different 
wedge angles. Very good agreement was found for the cases 
of Of- 15** and a =18°. Here, the pressure coefficient is 
defined as [pvvaii^C^^^) pooWiJ and the heat-transfer coef- 
ficient as 


c 




8T 

K ~ — sec0; 

by 


r . P P 

P«Woo[(e+ “ + — )oo- (e+ — )^.] 

P 2 p 


For Of =15° there is no separation, while for of=I8°, the 
separation and reattachment are very close to the experiment. 
The general features of the computed results for the a = 24° 
case have the correct trend but the extent of separation is 
slightly underpredicted. Nevertheless, the plateau pressure 
and the magnitude and location of the peaks of surface 
pressure and heat transfer are correctly predicted. Note that, 
in the experiment the peaks of pressure and heat transfer 
coincide, and in the computed results the peak of heat transfer 
is ahead but very close to those of surface pressure. This is 
also observed in the case of a= 18°. 

With the entire flowfield calculated, we can examine the 
details of the interaction between viscous and inviscid flows. 
Fig. 9 shows the isobar distributions for wedge angles of 15°, 
18°, and 24°. One of the most striking features is that the 
static pressure is neither constant across the boundary layer 
nor constant along the simple straight characteristic lines, as 
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Fijj. 10 a) pressure varialion and b) w-velocity boundary layer 
profiles near separation for « = 18^. 


treated by Myring.*^ Figure 10 shows the static pressure 
variation and the corresponding w-velocity profiles across the 
boundary layer at locations near separation for a = 18 . The 
pressure can be treated as constant only very near the wall, 
and then it decreases drastically to the edge of the boundary 
layer. The total variation may be as much as a factor of two or 
more, depending on the location and wedge angle. 
Presumably, the normal pressure gradient develops as a result 
of the curvature of the streamlines in the boundary layer. This 
curvature can be generated by the response of the boundary 
layer to the longitudinal pressure gradient or form the cur- 
vature of the surface upon which the boundary layer is 
growing. The pressure rise is first detected near the wall and, 
because of the large Mach number, takes quite a distance to 
reach the edge of the boundary layer. The ratio of this dis- 
tance to the boundary-layer thickness may be very large, and 
hence, (A/? /Ay) > {Ap /Ax) may occur. 

Downstream of the corner, the compression waves_coalesce 
into a shock wave. The intersection of the leading edge shock 
with this induced shock shows another interesting feature. 
Here both the leading edge shock and the induced shock are 
right-runing waves. The intersection is classified by Edney 


RESULTANT SHOCK , 




P/Po 


Fig. 12 Pressure profiles at three locations downstream of in- 
tersection point of leading edge and induced shocks. 


as Type VI interference in which, after the intersection of two 
shocks, a single stronger resultant shock, and expansion fan, 
and a slip surface are formed (Fig. 1 1). Along the streamline 
B, through multiple compression, the pressure in region 3 is 
higher than the pressure in region 5, which is compressed 
along streamline A by a single resultant shock. Behind the ex- 
pansion fan the pressure and flow direction in region 4 are the 
same as those in region 5, which is separated from region 4 by 
a slip surface. For a = 15% the expansion fan does not reach 
the wedge surface within the computational region. The wall 
pressure monotonically increases and the peak of maximum 
pressure is expected to lie downstream of the computational 
domain (Fig. 6a). For a = 18% the influence of the expansion 
fan on the wedge surface is observed in the computational 
domain (Fig. 9); the surface pressure and heat transfer reach 
their maximum and then decrease gradually (Fig. 7). Since the 
intersection angle is small, the expansion fan is so weak that it 
can barely be detected. As the wedge angle increases further to 
a -24% the correlation of the interference of leading edge 
shock and the peak pressure is evident (Fig. 9). (The dashed 
line represents, approximately, the trajectory of the leading 
edge shock.) The high pressure in region 3 (Fig. 11) is at- 
tributed to the viscous effect of smoothly bending the 
streamlines while passing the compression corner, and the 
high heat -transfer rate is attributed to the high temperature 
and the local thinning of the boundary layer behind the com- 
pression weaves and the induced shock. The expansion fan 
causes a rapid decrease in pressure. Figure 12 presents the 
pressure profiles at three different locations, AA , BB , CC 
(indicated in Fig. 9), and shows the details of the pressure 
field change from overcompression, through expansion, to a 
state of overexpansion relative to the inviscid solution. The 
expansion fan bounces back and forth, partly transmitted and 
partly reflected from the viscous layer and slip surface, bet- 
ween the wall and the external resultant shock, and the flow 
asymptotically relaxes to a state of normal weak interaction. 
It is interesting to notice that, for a flow at Mach number 14.1 
over a 24® wedge, the pressure rise,/? is 225 for the isen- 
tropic solution, and 58 for the inviscid Rankine-Hugoniot 
relation, whereas (pp^)^^^ = l02 for the viscous flow 
solution in the present calculation. 

Indeed there exists no discontinuity or slip surface in the 
viscous flow; the slip surface between regions 4 and 5 (Fig. 1 1) 
is smeared into a shear layer. Figure 13 shows the density con- 
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Fig. 13 Density contours and details of shear layer and boundary 
layer. 


Fig. 14 Velocity, tem- 
perature, and density profiles 
across boundary layer and 
shear layer. 
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Fig. 15 Locations of leading edge shock, induced shock, resultant 
shock, edge of boundary layer, boundary-layer displacement 
thickness, and line of zero u-velocity for a = 24“ . 


tours and clearly displays the trajectories of the leading edge 
shock, resultant shock, shear layer, boundary layer, and the 
cold-wall high-density sublayer. Figure 14 shows the detailed 
plots of w-velocity, density, and temperature across the boun- 
dary layer and the shear layer at the location {x/L) = 179 for 
a =24 . Total pressure loss is larger in region 5 than across 
the single resultant shock; hence, the temperature is higher, 
and the velocity and density are lower in region 5 than those in 
region 4. 

Figure 15 indicates the locations of leading edge shock, in- 
duced shock, resultant shock, edge of boundary layer, boun- 
dary-layer displacement thickness 6*, and the line of zero 
w-velocity for the case of a = 24”. Here, the edge of the 
boundary layer is defined as the smallest y for which 
A(pw) /pw<0.02. The locations of the peak heat-transfer rate 
and the peak wall pressure are labeled A, and A 2 , respec- 
tively. The boundary-layer displacement thickness 5* at first 
increases as pressure increases, up to the corner, then 
decreases as pressure continues to increase, reaching a neck 
region at about the point of peak pressure, and then starts to 
increase as the normal state of the boundary layer resumes. 
The thickness of the boundary layer before and after the com- 
pression are different by about one order of magnitude. 

IV. Conclusion 

An efficient time-splitting finite difference scheme has been 
used to obtain steady-state solutions of the Navier-Stokes 
equations for supersonic and hypersonic laminar flows over a 


compression corner. Favorable comparisons with previous 
calculation and with experiment indicate that the present 
calculations are accurate. The pressure profiles are neither 
constant across the boundary layer nor constant along simple, 
straight characteristic lines, as has been assumed in some 
analyses. Surface pressures higher than predicted by inviscid 
shock wave theory are obtained, due to the compression of the 
smooth bending of the streamlines. Edney Type VI in- 
terference of the leading edge shock with the induced shock 
results in an expansion fan which produces a large peak in the 
surface pressure and heat transfer. Consequently, it is im- 
portant to include the leading edge shock wave in the present 
study of such hypersonic flows. 
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